Lifetime of a massive particle in a de Sitter universe 
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We study particle decay in de Sitter space-time as given by first order perturbation theory in 
an interacting quantum field theory. We show that for fields with masses above a critical mass m c 
there is no such thing as particle stability, so that decays forbidden in flat space-time do occur there. 
The lifetime of such a particle also turns out to be independent of its velocity when that lifetime is 
comparable with de Sitter radius. Particles with lower mass are even stranger: The masses of their 
decay products must obey quantification rules, and their lifetime is zero. 
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Some important progress in the astronomical observa- 
tions of the last ten years [l|, l2| have led to the surpris- 
ing conclusion that the recent universe is dominated by a 
" dark" exotic form of energy density that acts repulsively 
at large scales. The simplest and best known candidate 
for the "dark energy" is the cosmological constant, and 
the de Sitter geometry, which is the homogeneous and 
isotropic solution of the cosmological Einstein equations 
in vacuo, appears to take the double role of reference ge- 
ometry of the universe, namely the geometry of spacetime 
deprived of its matter and radiation content and of the 
geometry that the universe will approach asymptotically. 

One might think that the presence of a cosmological 
constant, while having a huge impact on our understand- 
ing of the universe as a whole, would not influence micro- 
physics in its quantum aspects. This is also the viewpoint 
taken in the context of inflationary models 0] , where the 
effective cosmological constant is many orders of magni- 
tude larger than the one observed today. However this 
conclusion may have to be reassessed. Indeed, in pres- 
ence of a cosmological constant, however small, it is the 
notion of elementary particle itself which has to be re- 
considered, since the usual asymptotic theory is based 
on concepts which refer closely to Minkowski spacetime 
and to its Fourier representation, and do not apply to the 
de Sitter universe which is not asymptotically flat; in fact 
a true asymptotic theory does not exist at present for de 
Sitter space. A possible basic approach is perturbation 
theory; unfortunately, calculations of perturbative am- 
plitudes which in the Minkowskian case would be simple 
or even trivial become rapidly prohibitive or impossible 
in the de Sitter case: this in spite of the fact that one is 
dealing with a maximally symmetric manifold. 

In this letter we have tackled one such calculation, 
namely that of the mean lifetime of dS unstable scalar 
particles. The results exhibit significant differences com- 
pared to the Minkowski case and decay processes which 
are normally forbidden become possible (as exhibited by 
Nachtmann Q in a special case) and vice versa, processes 
that are normally possible are now forbidden. The maxi- 
mal symmetry of the dS universe allows for the introduc- 
tion of a global mass operator, one of the two Casimir 
operators of the dS group SO(l,d) (see e.g. ||); this 
quantity is conserved for dS invariant field theories and 



it still makes sense to follow Wigner Q in associating a 
particle with a unitary irreducible representation of the 
dS group labeled by a mass parameter, as we do here. 
However, in contrast with the Poincare group, the ten- 
sor product of two unitary irreducible representations of 
masses mi and m% decomposes into a direct integral of 
representations whose masses m do not satisfy the "sub- 
additivity condition" m > mi + mi', all representations 
of mass larger than a certain critical value (principal se- 
ries) appear in the decomposition. This fact was shown 
in [6| for the two-dimensional case and will be established 
here in general. This means that the de Sitter symmetry 
does not prevent a particle with mass in the principal 
series from decaying into e.g. pairs of heavier particles. 
This phenomenon also implies that there can be nothing 
like a mass gap in that range. This is a major obstruc- 
tion to attempts at constructing a de Sitter S-matrix; 
the Minkowskian asymptotic theory makes essential use 
of an isolated point in the spectrum of the mass oper- 
ator, and this will generally not occur in the de Sitter 
case. We will also show that the tensor product of two 
representations of sufficiently small mass below the crit- 
ical value (complementary series) contains an additional 
finite sum of discrete terms in the complementary series 
itself (at most one term in dimension four). This implies 
a form of particle stability, but the new phenomenon is 
that a particle of this kind cannot disintegrate unless the 
masses of the decay products have certain quantized val- 
ues. Stability for the same range of masses has also been 
recently found [|[ in a completely different context. 

We will resort to first order perturbation theory in our 
calculations. These arc made trivial in the Minkowski 
case by the use of momentum-space, but this is not so 
in the de Sitter space. We restrict at first our attention 
to the principal series and start by deriving a new gen- 
eral formula expressing the decay probability of a particle 
which applies to both Minkowski and de Sitter space- 
times in dimension d, where the de Sitter manifold is 
identified with the hyperboloid {x <E R d+1 : x 2 — 
Xq — x\ — . . . — x 2 d = — R 2 } in the (d + l)-dimensional 
Minkowski space. Differences will come in later. 

A Klein-Gordon neutral scalar field <f> with mass m > 
is characterized by its two-point vacuum expectation 
value w m (x, y) = (O, <j){x) 0(y)f2) which is uniquely 



2 



specified up to a constant factor by requiring invari- 
ance, locality, and a suitable spectral condition (having a 
thermodynamical physical interpretation in the dS case 
9, 10]). w m allows for the reconstruction of the Fock 
space of the theory and of a representation of the invari- 
ance group whose restriction to the one-particle subspace 
is irreducible and labeled by to. In the dS case to can be 
related to a dimensionless parameter v as follows 
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The range to > m c = [d—\)/2R (i.e. v real) corresponds 
to the principal series while < to < m c corresponds to 
the complementary series (y imaginary). These restric- 
tions ensure that w m is positive definite and therefore a 
quantum theoretical interpretation is available. Consider 
now an interaction 



jg{x) C(x) dx, £(x) = : M^M*)" 1 ■ ■ ■ <t> N {x) qN 



between 1 + N independent scalar fields </>o > </>i j • ■ ■ i 4>n 
with masses too, toi, . . . , to^t; self-interactions C(x) = : 
4>o(x) n : are a special case of this coupling. The spacetime 
dependent " switching-on factor" g(x) is there to take care 
of the infrared divergence of the integrals and amounts to 
putting the system in a box and allowing for a finite time 
duration of the interaction. In the end however g should 
be made to tend to 1 everywhere (adiabatic limit). 

Let \P = j dx f(x) (fio(x)Cl be a (normalized) one- 
particle state created by 4>o from the vacuum; f(x) en- 
codes the physical details about the quantum state of the 
unstable particle whose disintegration we aim to study. 
At first order in perturbation theory Wick's theorem 
gives the transition probability r lo;gi qN from ^ to any 
possible state containing qi + q 2 + . . . + qn particles cre- 
ated respectively by the fields <f>\, <f>2, ■ ■ ■ 4>n'- 



T = 7 J dxdudvdyf(x) f(y)g(u)g(v)K(x,u,v,y),(2) 

N 

K(x,u,v,y) = w mo (x,u) q \ w mj (u, v) q] w mo (v,y). 

3 = 1 

The four-point kernel K admits an easy graphical inter- 
pretation in a;-space as a bubble with two external legs. 
The product of two-point functions (the bubble) may be 
replaced by its K alien- Lehmann type representation: 

N 

Yl w mj (u, v)* = / da 2 p(a 2 \ to 1; ..., m n )w a (u,v). 

.(3) 

It is far from obvious that such representation exists; 
besides, calculating p is easy only in the Minkowski case 
for n — 2 (see below); difficult or impossible in other 
cases. While the amplitude @ is infrared divergent in 
the limit g — > 1, it is possible to replace just one of the g's 
by 1 in the integral. Then one integration can be done 



by means of the projector identity (dy denotes either the 
Poincare or the dS invariant measure) : 



w m {x, y) w m > (y, z) dy = C(m)5(m 2 - m' 2 )w m (x 1 z) 



Here the global structure of either the Minkowski or 
the dS spacetime enters crucially. For dS theories this 
identity holds only for the principal series; C(m) = 2tt 
(Minkowski); C{y) = 2n\ coth(7w)| (dS). There follows a 
general formula for the transition probability: 
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( x , y)f(y)dxdy 

p(ml;mi, . . . ,m N ). 



(4) 



Here F(x) = J w ma (x, y)f(y)dy is the wavefunction 
associated to /; the denominator is the squared norm 
of no longer assumed to be one. This formula has 
an interesting simple structure: the first factor does not 
depend on the number or nature of the decay particles 
but only on the wavefunction of the incoming unstable 
particle. The infrared problem is contained in this factor 
and has to be overcome when letting the remaining g{x) 
tend to 1. The second factor is just the relevant Kallen- 
Lehmann weight times the right combinatorial factor. 

We now focus on the decay of a particle of mass Too 
into two particles of mass Toi and briefly discuss the well- 
known Minkowski case first. The weight /)(ra§;mi,TOi) 
can be computed by Fourier transforming w 2 ni (x, y): 
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The appearance of Heaviside's 8 function forbids the de- 
cay of a particle into two that are globally heavier. This is 
a familiar consequence of the Poincare invariance of the 
theory. As for the adiabatic limit, the common choice 
is to choose g{x) as the characteristic function of some 
time interval T. It is then found that the transition prob- 
ability ([?]) is proportional to T and thus diverges when 
T — > oo. Fermi's golden rule tells us that the transition 
probability per unit time (see e.g. (llj ) has a finite limit 
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(27r) 7 2 JW)- 1 \f(p)\ 2 S(p 2 - m 2 )6(p°)dp 
J\f(p)\ 2 S(p 2 -m 2 )e(p°)dp 



(6) 



f(p) is the Fourier transform of the wavepacket f(x). 
The crucial factor (2p°)~ 1 in the numerator controls the 
dependence of the result ^ on the wavepacket /. In 
particular, the decay rate of a particle at rest in our frame 
can be obtained by letting |/(p)| 2 tend to S(p). In this 
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limit the factor (2p°) 1 becomes (2mo) 1 and one gets: 

#(m 2 -4m 2 ). (7) 
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Had we chosen the wavepacket of a particle with sharp 
momentum p, we would have obtained an extra Lorentz 
factor and the lifetime would be longer: 



t(v) = tq(1 — v 2 /(?) 2 , v — cp/p° 



(8) 



In the de Sitter case we use the dimensionless parameter 
v (see Eq. [IJ to label the two-point functions; they are 
proportional to Legendre functions of the first kind: 
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x (c 2 -i)-^p:S„(0; (9) 



z, z' belong to suitable tubular domains of the complex 
de Sitter spacetime; £ — z ■ z 1 / R 2 in the ambient space- 
time sense (see @ for details) . For d = 2 the positivity 
of piv 2 ; v, . . . , v) has already been established [4j for the 
disintegration into an odd number of particles having the 
same mass as the unstable particle. Here we need to ac- 
tually compute the Kallen-Lehmann weight i.e. to obtain 
p(k 2 ; v, v) = Pv{k) such that 
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Due to (|9|) , the generalized Mehler-Fock theorem [12j tells 
us that p v {n) can be expressed as the integral 
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which is well defined for masses such that | \mv\ < ^p-; 
this includes the principal series and a portion of the com- 
plementary series. Computing (llip involves only trigono- 
metric functions for odd d [13j. Remarkably, Mellin 
transform techniques [3] allow the computation for any 
dimension d (details will be published elsewhere): 
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Contrary to ([5]), the weight /? never vanishes. This means 
that for m > m c decay processes into heavier particles 
are always possible and thus, in that range of masses, 
one is not allowed to draw conclusions about the stabil- 
ity of a certain particle just from its being the lightest 



in a hierarchy. The Minkowskian result ([5]) is however 
recovered in the limit of zero curvature (R — > oo) that 
is achieved by setting k — moR and v = m\R. Lowest 
order corrections to the flat case give: 
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Am = mo — 2m\. The lack of particle stability (Am < 
0) is exponentially small in R. If Am > there is a 
correction to the flat case of the order of the cosmological 
constant A = ^ d ~^a~ 2 ^ ■ In the four dimensional case 
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All these effects are of course extremely small with 
the current value of the cosmological constant. What 
about particle physics at inflation? At that epoch mR ~ 
m x 10~ 15 GeV <C § for every particle of reason- 
able mass. Our results should therefore be extended 
to the remaining portion of the complementary series 
^-p < | Im^| < where all scalar particles lie at the 
inflation era (but: there is no complementary series in 
the Fermionic case). By analytic continuation of (fT2|) in 
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where p = (d — l)/2. The number of discrete terms is 
the largest N satisfying N < 1 + | Im v\ — p/2, or if this 
is negative. A particle of the complementary series with 
parameter k = i/3 can only decay into two particles with 
parameter z/=|(|/3|+/i + 2n), where n is any integer 
such that < 2n < p — \(3\, and the decay is instanta- 
neous. A particle with mass m <C m c can only decay into 
two particles of mass mi ~ m/V2. Even if the geometry 
of the universe at inflation was not exactly de Sitterian, 
this example indicates that quantum field theoretical ar- 
guments concerning particle physics at inflation might 
need revision. 

In trying to interpret the above results one can won- 
der whether they might be due to the thermodynamical 
properties [l(| of the fundamental state we have been 
using. We have tested this possibility against a similar 
computation in flat thermal field theory that however 
does not exhibit this phenomenon in two-particle decays. 
Another issue has to do with energy conservation and 
the relation mass/energy. dS invariant field theories ad- 
mit ten conserved quantities (in d — 4). The identifi- 
cation of a conserved energy among these quantities has 
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proven to be useful in classical field theory [15[ . The same 
quantity remains exactly conserved also at the quantum 
level although it becomes an operator whose spectrum is 
not positive Q even when restricted to the region where 
the corresponding classical expression is positive [1 5| ; the 
thermodynamical properties of dS fields arise precisely in 
this restriction Q. Energy is conserved also in the de- 
cay processes that violate mass subadditivity, once the 
adiabatic limit has been performed. The breakdown of 
the subadditivity property of masses in dS spacetime just 
reflects the nonexistence of an Abelian translation group 
and thereby of a linear energy-momentum space. 

We now consider the adiabatic limit problem and its 
meaning in the de Sitter context, in the case when all 
particles are in the principal series. A first complication 
is the existence of several choices of cosmic time, hav- 
ing different physical implications and the result might 
depend on one's preferred choice. In the closed model, 
the cosmic time t is related to the ambient space coor- 
dinates as follows: x° — R sinh (t/R). In strict anal- 
ogy to the Minkowski case, g(x) can be chosen as the 
indicator function of some cosmic time interval T, say 
g(x)=g T (x)=6(T/2-\t\). 

In the flat model the situation is a bit more tricky. 
Cosmic time is now defined by the relation x a + x d = 
R exp(i/i?); flat coordinates cover only half of the de 
Sitter manifold, namely all the events such that x° + x d > 
0. If we introduce the characteristic function hx{x) — 
9{Re T / 2R -x°- x d ) 6(x° + x d - Rer T / 2R ) then we have 
to add the contribution coming from the other half, i.e. 
g{x) — gr{x) — hr{x) + hx{—x). With these premises 
we have found that in both models the first factor in (|4|) 
diverges like T; thus it has to be divided by T to extract 
a finite result which is the same in both models: 

Hm ^Cjn) Jg(x)\F(x)\ 2 dx = 7 2 7rcoth(^) 2 
T ^°°Tff(x)w K (x,y)f(y)dxdy M 

Here the second (unforeseen) result comes in: in con- 



trast to the Minkowskian case the limiting probability 
per unit of time does not depend on the wavepacket! 
This result seems to contradict what we see everyday 
in laboratory experiments, a well known effect of special 
relativity (Eq. [7]). Furthermore, in contrast with the vi- 
olation of particle stability that is exponentially small in 
the de Sitter radius, this phenomenon does not depend 
on how small is the cosmological constant. How can we 
solve this paradox and reconcile the result with every- 
day experience? The point is that the idea of probability 
per unit time (Fermi's golden rule) has no scale-invariant 
meaning in de Sitter: if we use the limiting probability 
to evaluate amplitudes of processes that take place in a 
short time we get a grossly wrong result. This is in strong 
disagreement with what happens in the Minkowski case 
where the limiting probability is attained almost imme- 
diately (i.e. already for finite T). Therefore to describe 
what we are really doing in a laboratory we should not 
take the limit T — > oo and rather use the probability per 
unit of time relative to a laboratory consistent scale of 
time. In that case we will recover all the standard wis- 
dom even in presence of a cosmological constant. But, if 
an unstable particle lives a very long time (>> R) and we 
can accumulate observations then a nonvanishing cosmo- 
logical constant would radically modify the Minkowski 
result and de Sitter invariant result will emerge. This 
result should not be shocking: after all erasing any in- 
homogeneity is precisely what the quasi de Sitter phase 
is supposed to do at the epoch of inflation; in the same 
way, from the viewpoint of an accelerating universe all 
the long-lived particles look as if they were at rest and so 
their lifetime would not depend on their peculiar motion. 



We thank T. Damour, M. Gaudin, G. Gibbons, M. Mil- 
gram and V. Pasquier for enlightening discussions. U. M. 
thanks the SPhT and the IHES for hospitality and sup- 
port. 



[1] A. G. Riess et al., Astronomical Journal 116, 1009 (1998). 

[2] S. Perlmutter et al., Astrophys.J.517:565 (1999) 

[3] A. Guth, Phys. Rev. D23, 347 (1981). A. Linde, "Particle 

Physics And Inflationary Cosmology," Chur: Harwood 

(1990). 

[4] O. Nachtmann: Sitz. Ber.Ost. Ak. d. Wiss. II 176, 363 
(1968). 

[5] F. Giirsey, Introduction to the de Sitter group. In: Group 
theoretical concepts and methods in elementary parti- 
cle physics pp. 365-389 Gordon and Breach, New York 
(1964) 

[6] Y. A. Verdiyev, "Harmonic Analysis on Homogeneous 

Spaces", Hadronic Press, (1986). 
[7] E. Wigner, Ann. of Math. (2) 40 (1939), no. 1, 149-204. 
[8] K. Skenderis and P. K. Townsend, J. Phys. A 40, 6733 

(2007). 

[9] J. Bros, J. P. Gazeau and U. Moschella, Phys. Rev. Lett. 



73 (1994) 1746. J. Bros and U. Moschella, Rev. Math. 
Phys. 8 (1996) 327. J. Bros, H. Epstein and U. Moschella, 
Commun. Math. Phys. 196, 535 (1998) 
[10] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15 
(1977) 2738. 

[11] M. Veltman, "Diagrammatica" , Cambridge University 
Press, (1994). 

[12] W. Magnus, F. Oberhettinger, and R. P. Soni, "Formulas 

and Theorems for the Special Functions of Mathematical 

Physics", Berlin: Springer- Verlag (1966). 
[13] M. Gaudin and V. Pasquier, unpublished 
[14] O. I. Marichev: "Handbook of Integral Fransforms 

of Higher Transcendental Functions". Chichester: Hor- 

wood, 1982 

[15] L. F. Abbott and S. Deser, Nucl. Phys. B 195, 76 (1982). 



